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Abstract 

We reexamine the one-loop corrections to the mass of the pseudo-scalar Higgs 
boson, using the effective potential. In the absence of the chargino and neutralino 
contributions its mass exhibits a large scale dependence in the large M^/ 2 regime, 
especially for values of tan (3 > 20. Thus, although of electroweak origin, the 
heaviness of the M1/2, in conjunction with the largeness of tan (3, makes these 
corrections very important for establishing a scale independent result and an un- 
ambiguous determination of the pseudo-scalar mass in this region of the parameter 
space. 



1 Introduction 



The Higgs sector of the MSSM has been put under experimental scrutiny the last few 
years and reliable lower bounds on the Higgs masses have been already established im- 
posing tight constraints on supersymmetric models 0. 

Radiative corrections to the Higgs masses have been extensively studied, especially 
for the lightest CP-even Higgs state for which a theoretical upper bound has been 
established of ~ 130 GeV 0-15]- Among the others states the CP-odd Higgs A may play 



a significant role in processes in which CP is violated pq] , and also in cosmology since 
the lightest supersymmetric particle (LSP) pair annihilation to a fermion pair through 
A exchange may enhance the corresponding cross section for large tan/3, yielding LSP 



relic densities compatible with the recent astrophysical data [17]. The latter process is 
sensitive to the mass of the pseudo-scalar Higgs and a reliable determination of its mass 
is highly demanded. 

Existing studies in literature discuss the radiative corrections to the pseudo-scalar 
Higgs mass using either the effective potential approach or by studying the loop correc- 
tions to the propagator |2|^.[To Hl^1 . The mass in the second approach is the physical, or 
pole, mass rriA.{pole) which is certainly scale independent and differs from that calculated 
using the effective potential denoted by tua- The two masses are related by |||| 

m 2 A (pole) =m 2 A + Il AA (0) - U AA (m A (pole)) , (1) 

where Uaa{p 2 ) are the corrections to the pseudo-scalar propagator at momentum p. In 
some studies it is claimed that rriA is scale independent, however this assertion is not 
entirely correct since the difference n^^m^poZe)) — n^^O) includes small logarithmic 
parts, which depend explicitly on the scale Q. For completeness these should be added 
to m A to render scale invariance. 

In this note we consider all radiative corrections to the pseudo-scalar mass as deriva- 
tives of the one-loop effective potential. We particularly show the importance of the 
chargino and neutralino corrections in establishing a result which is scale independent 
and approximates accurately the pole mass. Long known third generation sfermion con- 
tributions are not by themselves adequate to yield a scale independent result for large 
Mi/2 values. The reason is that the Renormalization Group Equation (RGE) of the 
Higgs mixing parameter m|, entering into the expression for ttla, includes terms 

' ii (-6g 2 2 M 2 - , (2) 
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resulting to a scale dependence of m 2 not cancelled by the third generation sfermion 
corrections alone. In the large M\n and large tan/3 regime such corrections are not 
numerically small. 
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In our approach for the determination of the pseudo-scalar mass we have duly taken 
into account all contributions, including the chargino/neutralino corrections, and have 
observed that the latter contribute significantly to the stabilization of the pseudo-scalar 
Higgs boson mass with respect the scale Q. If these contributions were neglected stabi- 
lization would be spoiled in the large M1/2 region, where the running of the parameter 
m\(Q) due to the gauginos becomes important. The situation would be even more dra- 
matic if in addition to having large My 2 we are in the large tan/? regime where such 
deviations from stability are enhanced as being proportional to tan (5. 

From Eq. (|]) we see that the effective potential mass differs from the physical mass 
by 11(0) — U(m\). In our scheme, and for a more reliable estimate of the pseudo- 
scalar mass, besides the chargino and neutralino corrections we have also included the 
contributions of the remaining sectors to rriA as well as the leading logarithmic corrections 
of 11(0) — H(m\). Including these the resulting mass is certainly scale independent, to 
this loop order, and is found to be very close to the one-loop pole mass. Differences 
are less than 2%. In this way we have approximated satisfactorily the pole mass and 
avoided the complexities of calculating one-loop integrals, which are usually expressed 
by the Passarino-Veltman functions. 

2 Radiative Corrections to the pseudo-scalar Higgs 
Mass 

We assume a low energy supersymmetric theory described by the superpotential 

W = h t Hje QU C + h b H?e QD C + h T Hfe LE C + fiH^e H x , (3) 

where the elements of the antisymmetric 2x2 matrix e are given by €\2 = —£21 — — !■ in 
it only the dominant Yukawa terms of the third generation are assumed nonvanishing. 
The soft-SUSY breaking part of the Lagrangian is given by, in an obvious notation, 

Acalar = ~ ^m^|0;| 2 - {m\H H 1 + fl.c) 

i 

- (A t h t H^e QU C + A b h b Hfe QD C + A T h T H?e LE C + h.c) 

- - (M X BB + M2W (i) W {i) + M 3 GG + h.c.) . 

At the tree level the scalar potential of the theory is given by 

V° = mllH^ + mllH^ + mKH^H^ + h.c) 

+ (^^)(I^ I 2 -I^ 2 T) 2 , (4) 

where = m 2 H . + [i 2 . In this equation we have only kept the part that depends 
exclusively on the neutral components of the Higgs fields H® 2 which is relevant for the 
electroweak symmetry breaking. 
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The one-loop correction to the effective potential is 

f 



^ 1 ^E(-^'(2. SJ + i )m j(.„9l-|) 



(5) 



where mj are field dependent masses and sj we denotes the spin of the J-particle. 

The minimization of the one-loop corrected effective potential V 1 = V° + AV 1 yields 
the following relations 

2 8 raj — ml tan 2 j3 



sin 2(3 = 
where we have defined [] 



2m 2 

— 2 i — 2 ' 



# 2 + #' 2 tan 2 (3 - 1 



m,- = m,- + S 



a^ 1 



' ' Si " <9(Re#°) 2 



(6) 



(7) 



As discussed in the introduction the relation between the pole and the mass calculated 
using the effective potential is given by Eq. ([J). In this m\ is the non-zero eigenvalue 

where </>,■ = Inxf/p. At the tree level 

<*?> 



of the matrix ■= 



1 d 2 V° 



1 <9 2 ^° 

2 d6id<b. 



(ml + ffin) cos /3 sin /3 



tan /3 1 
1 cot/? 



(8) 



The matrix above has one zero eigenvalue corresponding to the mass of the Goldstone 
mode which is eaten up by the Z-boson. The other eigenvalue is the mass squared of 
the pseudo-scalar Higgs, given by m\ = m\ + m\ in the lowest order. 
At the one-loop the Eq. (0) is modified to 



i d^v 1 



2 d(f>id(j)j 



T,2 



™2> 



m 



cos /3 sin /3 



tan /3 
1 



1 

cot /3 



+ 



x 



647T 2 ^ 



-I) 2 - (2 Sj + 1) m 2 (ln^|-l 



d 2 m 2 



d 2 



d 2 m 2 d 2 m 2 
d4>2d<t>i dff) 2 . 



J 



0<p\O<p2 

- 2 



dm 2 



(9) 



d(ReH$) 2 



(«y> 

where we have used the one-loop minimization conditions of Eq. @ as well as the 
relations of Eq. (|7|) and the one-loop corrections to the scalar potential given by the 
Eq. (|5|). This is the master formula we are going to use for the calculation of the one- 
loop corrections to the CP-odd Higgs boson in the effective potential approach. 

Using the minimization condition relating (m\ + mf) to the Higgs mixing parameter 



(see Eq. 



the above formula can always be cast into the form 

1 d 2 V l 

2 dfadfa / H o\ 



-(ml + A) 



tan (3 1 
1 cot/? 



(10) 



^ote that in our notation v, = v 1 = cos/3, v 2 = -75 sin/?, = g 2 (v 2 + v\)j2 = g 2 v 2 /A 



3 



resulting to a pseudo-scalar mass squared m A given by 



_ 2(m§ + A) 
mA ~ sin 2(3 ■ 

For the calculation of the CP-odd Higgs mass we thus need calculate the contribution of 
each particle species to A. The dominant third generation sfermion contributions to this 
quantity have been long known. The neutralino and chargino contributions are claimed 
to be small and hence not expected to yield substantial corrections to the pseudo-scalar 
mass [fLU[| . However this may not be true when the soft gaugino mass M 1( / 2 is large as 
pointed out in the introduction. In this regime such corrections may be sizeable and 
should be duly taken into account. 

It is well known that radiative corrections to rriA, calculated through the effective 



potential, are not stable with changing the renormalization scale (see for instance Ref. |T9 



and references therein). Empirically, one calculates m A at an average stop mass scale, 
Qt ~ \J m h m t 2 i m which case the radiative corrections of the third generation sfemions 
are small and can be safely neglected. In this case the pseudo-scalar mass squared is given 
by m A = —2m\{Qi)/sm2f3{Qi), where m\ is the Higgs mixing parameter. Although in 
principle this is correct, the contributions of charginos and neutralinos are not small 
at this scale, especially when M1/2 is large. In some cases this may produce an error 
in the determination of its mass as large as 25%. Excluding the chargino/neutralino 
contribution is legitimate provided the relevant scale is not taken to be Q^, but rather 
the average of the chargino and neutralino masses defined by the following expression 

Q\ = \ ( (m%) + (ml) ) . (12) 

In the definition above (m 2 ^) , ( m \) denote the averages of the chargino and neutralino 
masses squared. At this scale it is legitimate to neglect their contributions. This scale 
however may differ substantially from when M1/2 is large. Therefore one expects large 
variations of the between these two scales, if only the third generation sfermion 

contribution is kept, resulting to a poor determination of the pseudo-scalar mass. In the 
following we shall calculate all corrections to the pseudo-scalar mass. The gauge and 
Higgs boson contributions to A although less important should be also included to yield 
a result that is scale independent and approximates satisfactorily well the pole mass. 
An alternative way to calculate the quantity A is through the one-loop corrections 



to the pseudo-scalar propagator n^. In fact following Ref. [0 we find that A can be 
expressed as 

A = S/3 Cf3 ( U AA (0) - s 2 p - - c}^) (13) 

where sp = sin (3 , cp = cos (3 , and t\ y2 are the one-loop contributions to the v.e.v's 
(ifi i2 ) . The latter can be calculated by either using the effective potential [|l8j or by 
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employing diagrammatic techniques. In principle the two approaches yield identical 
results. The gauge dependence of the tadpole graphs cancels against gauge dependence 
of loop diagrams in Haa which involve gauge and Goldstone bosons. Lacking the form 
of the effective potential in gauges other than the Landau gauge, such as the popular t' 
Hooft gauge for instance, we prefer to rely on diagrammatic techniques to calculate the 
gauge and Higgs contributions to the tadpoles whenever needed. 

In the following we shall discuss the various contributions to the quantity A. To 
establish our notation we briefly recall the results of Refs. 0, U , and we first consider 
the third generation sfermion contributions to the CP-odd Higgs scalar. The one-loop 
corrections are given by Eq. (|10D with A given by 

1 



A q 



327T 2 



N f h}A 



f 



f=t,b,r h 



m 



h 



[/(r4)-/(m|)]. 



In the equation above Nf is the color factor and the function f(m 2 ) is defined by 

,2 



f{m 2 ) = 2m 2 



i m i 

V" 1 



(14) 



The remaining sfermion contributions are zero because of the vanishing of the corre- 
sponding Yukawa couplings. 

For the chargino contribution we need their field dependent mass matrix A4c • Its 
squared M? c — A4 c Aic is given by the following relation 



M 2 + g 2 \H° 2 \ 2 -g(M 2 H<? + f,H$ 



-g{M 2 Hl* + ^H°) 



012 



From this, after a straightforward calculation, it follows that the chargino contribution 
to the quantity A is given by 

c _ g 2 M 2 fi 



A 



9 

16n 2 m% 



[f(m 2 di )-f(m 2 d2 )] 



(15) 



The contribution of the neutralinos to A is less trivial to be found. Their field 
dependent mass matrix can be put in the following form 



M 



N 





V2 
V2 





M 2 

V2 
V2 









\ 



















/ 



(16) 



whose squared is defined by 



M 



Its field dependent eigenvalues m% 
eigenvalue equation 



N 

a - 



M\M N 



(17) 



1, 2, 3, 4 are real and are determined from the 



h(X) = det(M 2 N - XI) = 0. 



(18) 
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Obviously m% are the physical neutralino masses squared when the Higgs fields are on 
their v.e.v's. The function h(X) can be written as 

h(X) = A 4 + AX 3 + BX 2 + CX + D , (19) 



where the coefficients A, B, C and D depend on the neutral Higgs fields Hf . 

In order to bring the second derivatives of the potential into the form of Eq. (|10|) 

dm 2 - dm 2 - 

we need the derivatives am and », f° where 6a is ImH?. The relevant derivatives 



d(RcH?) 2 



entering the Eq. (pi) are 



dm% 

Z a 



d(ReH, 



0\2 



h'(m% ) 



[ml/A, + (m%J 2 Bi + (mlJCi + D< 



2a' 



(20) 



and also 



d 2 m 2 - 



fc'K* ) 



2 \3j// 



(21) 



In these equations 



A = 



OA 



and the same holds for B, C and D. In Eqs. 

dh(X) 



and A'?. 



<9 2 A 



(22) 



h'(m% ) stands for the derivative 



h'(ml) 



dX 



(23) 



In these all Higgs fields are meant on their v.e.v's and hence all masses appearing are 
the physical tree level neutralino masses. In deriving Eq. (|21|) we have used the fact that 



A' = B' = C' = D' = 0, when 



0. 



(24) 



Using these one can find that the neutralino contribution to A is given by 



A z = 



327T 2 



a=l 



where 



F ( m l ) = TiJ^T 
Za h'imi 



Z a ' 



(m|J 2 c x + (m| ) c 2 + fx 2 M 1 M 2 c 3 f{m\ ) . 



z a 



Z a ' 



(25) 



(26) 



The analytic expressions for the derivative h'(m 2 ~ ) appearing in the equation above is 

Z a 

easily expressed in terms of the physical masses and is given by 



fc/(m|J = jQ (m| a — m| ) a, b = 1,...,4 



(27) 
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The quantities 01,2,3, are given by 

Cl = M 2 + tan 2 6 W Mi 

c 2 = -tan 2 6 W Mi (/i 2 + M|) - M 2 (/i 2 + M x 2 ) 

c 3 = Mi + tan 2 6 W M 2 . (28) 

For a complete analysis the contributions of the gauge and Higgs bosons to the quan- 
tity A, although small, should be also included. Their contributions can be evaluated 
more easily using Eq. flilf). In this case the needed corrections to the pseudo-scalar 
propagators and tadpoles can be read from Ref. [Q. The pseudo-scalar propagator at 
zero momentum transfer is given by 



2 , e 2„ A 



16 vr 2 IWO) = 9 - { 2 D(H + , W) + £>(#, + £>(/i, 
8 1 <r c z 



M 2 



c 2 



4, ( c 2 ^ £(Z, if) + 5», £>(Z, /») ) ] 



+ ^ [ c 2/3 c 2a ( /(if) - f(h) ) -3c% f(A) + (1 - 3 ^ ) f(Z) } 

+ \ 2 [ (s 2 4p - J ( 1 + 4p )) f(W) - c% /(if+) ] 
- 2 M 2 , D(H + ) W) - 8 - i f{Z) }. (29) 

In this expression 

s 2a = sin(2a) , c 2a = cos(2a) , s 2l3 = sin (2/3) , c 2/3 = cos(2/3) 

s Qj a = sin(a — ft) , c a/ 3 = cos(a — ft) , = sin(a + ft) , c a/ g = cos(a + ft) . 

The functions appearing in Eq. (p9|) are given by 

f(i) _ f( ? -) _ m 2 /(i) - m 2 /(,) 

D(i,j) = ^ , D(i,j) = L^dl (30) 

mf — mi mf — mi 

where f(i) = f(mf). The function / was defined in Eq. (|14j) . H, h, A, H + denote the 
heavy and the light CP-even Higgs, the pseudo-scalar Higgs and charged Higgs boson 
respectively. W, Z stand for the charged and neutral gauge bosons respectively. 

The combination of the tadpoles needed to calculate the contributions of the gauge 
and the Higgs bosons to the quantity A, through Eq. fll3"D , is given by Ref. |~2 



ti 2 h 2 

af^{-!^-<5 + S>^ + > + 

1 



+ ^? [ s l ( c a ~ 3 s l ~ s 2« tan/? ) + c| ( s 2 - 3 c 2 - s 2q cot/? )] f(h) 
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+ ^ [ 4 ( s « - 3 c « + tan/3 ) + c| ( c 2 - 3 s 2 a + s 2a cot/3 )] /(F) 



+ (#-^J /(Z) )' (31) 

With this we exhaust the list of all nonvanishing contributions to the quantity A. 

For a scale independent result we should also include the leading logarithmic contri- 
butions of the difference Haa (0) ~^-aa{'^a) involving the renormalization scale Q. The 
log(Q 2 ) pieces of this expression, which relates the pole and effective potential masses, 
are listed below. These can be derived from the analytic expressions found in Ref. [|1^ . 
We write the log(Q 2 ) dependent part of this difference as 



n Lo9 = [ IWO) - U AA (m A ) 



logQ 2 

- (£*)^lo g (f). (32) 

i 

The remaining part of IIaa(O) — 11^(^.4) does not involve logarithms of the scale Q 2 . 
The nonvanishing contributions to the quantities Xi appearing in the expression above 
are given below 

X fermions = 3 Cg h t + 3 S p /ij, 4" Sp fl T , 

y „2 Y _ 9 2 

char ginos y ; -^neutralinos 



cos 2 9w ' 
3 

Xgauge+Higgs = ~ 9 ( ^ + ^an /3 ) . (33) 

Within the logarithm in the Eq. (|32|) the pseudo-scalar mass appears. Any other mass 
can also lead to a scale independent result. However we have numerically verified that 
this choice is the most appropriate in the sense of better approximating the pole mass. 

3 Results and Discussion 

Following the previous discussion we see that for defining a scale independent pseudo- 
scalar mass we have to incorporate to the effective potential mass of Eq. (|ll|) the 



logarithmic terms given by Eq. (|32|) . Therefore we define mr A as 



rh A = m A + Yi Log . (34) 

This is scale independent at the one-loop order and its difference from the pole mass 
is expected to be small. We have scanned the MSSM parameter space assuming uni- 
versal boundary conditions at the Unification scale and have indeed verified that fh A 
is very close to the pole mass. Differences are less than 2%. We have also seen that 
the logarithmic corrections T\.L og are small and therefore not very significant. Thus the 
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bulk of the radiative corrections to txia is carried by the effective potential mass txia- A 
sample result depicting very clearly the situation is displayed in the Figure [TJ. For the 
given inputs we plot the mass friA as function of the scale Q (dashed line) and compare 
it to the one- loop pole mass (dashed-dotted line), in a is indeed very close to the pole 
mass and almost independent of the scale Q for values of Q > Q%- In the displayed 
figures the solid line corresponds to the value of the effective potential mass rriA where 
the contribution of charginos and neutralinos is omitted. For defmiteness we shall denote 
this by 77^(0). This is scale dependent and deviates from the pole mass. This deviation 
becomes significant for large values of Mi/ 2 ^> m Q and large tan/3. In some cases the 
deviation of myi(O) from the pole mass can be as large as 25% depending on the scale it 
is calculated. This shows the significance of the chargino and neutralino corrections in 
estimating the pseudo-scalar Higgs boson mass. Only at Q « this coincides with the 
pole mass. 

Figure 2 shows the same situation for the same values of the soft masses and trilinear 
couplings but for a smaller value of tan (3. The effect is less dramatic and the inclusion 
of the charginos and neutralinos, denoted by x m the sequel, amounts to no more than 
~ 3%. In Figure 3 we have taken Mi/ 2 to be comparable to mo but the angle tan/3 
is taken to be large (= 30). In this case the resulting x corrections are of the order of 
~ 7%. In all figures the mass m^(0), unlike friA, is scale dependent resulting to a poor 
determination of the pseudo-scalar mass. From the discussion above it becomes obvious 
that the effect is more pronounced in regions of the parameter space for which both M\/2 
and tan /3 are large. The reason is that in the large tan j3 regime m 2 A is given by 

m\ « -tan ,3(7713 + A). (35) 

Therefore the gaugino soft terms of Eq. (0) that enter the RGE of the parameter 
m|, which become important for large values of the parameter M1/2, give a large scale 
dependence if not cancelled by corresponding terms in the corrections A. Failing to 
include the x contributions results to incomplete cancellations of log(Q 2 ) terms in (|35p, 
and this effect is augmented for large values of tan/3. This is the reason the mass m^O) 
shows this drastic scale dependence. 

The masses friA and m^O) are almost equal at scales comparable to Qy. and thus 
both very close to the pole mass. This is expected since at Q ~ the contribution of 
charginos and neutralinos are small and the H-iog terms, as said before, do not contribute 
significantly. The smallness of the x's contributions at Q ~ is due to the fact that 
their masses are close to resulting to small logarithmic corrections within the expres- 
sion of the effective potential. Therefore neglecting the contribution of the charginos and 
neutralinos is legitimate provided that the relevant scale of the calculation is taken to 
be Qx rather than Qi. These two scales however may lie quite apart. In the displayed 
figures the equality of the two masses is at the point where the dashed and the dashed 
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curve cross. One sees that this scale is quite close to whose value is shown in the 
figures. 

Below Qy. and Qj, and as we approach Mz, the mass in a shows an unstable behaviour 
with changing the scale Q which is however milder than that of m^O). This was expected 
since for such values of Q, which lie lower than the t and x masses, large logarithms enter 
the one-loop corrections of the effective potential. At such scales Q only the inclusion of 
the two-loop corrections can ameliorate the situation and yield a stable result. 

We conclude that for an unambiguous determination of the pseudo-scalar mass, as 
this is calculated from the effective potential, the contribution of the charginos and 
neutralinos should be taken into account. Their contributions is significant in regions of 
the parameter space for which M1/2 and tan/3 are large. Neglecting their contributions 
can be legitimate only if the calculation is performed at a scale close to the average of the 
chargino and neutralino masses. Adding the small contributions of the remaining sectors, 
gauge and Higgs bosons, as well as, the small leading logarithmic corrections of the 
difference 11(0) — U(m A ), results to a mass which is scale independent and approximates 
accurately the pole mass. 
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Figure 1: The pseudo-scalar mass as function of the scale Q for the inputs 
displayed in the figure. Qx{Qi) is the average chargino/neutralino (stop) mass. 
The solid line is the effective potential mass where only the third generation 
of sfermion contributes. In the dashed line the contribution of all species is 
taken into account as well as the the leading log wave function renormalization 
contributions. The dashed-dotted line is the pole mass. 
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Figure 2: The same as in Figure ffl with tan/? = 10. 
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Figure 3: The same an in Figure [l] with M1/2 to be comparable to mo- 
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